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OA . It is by now well known that the Boltzmann-Gibbs (BG) entropy Sbc ~ — X^J^li Pi Inpi can 

^ ' be usefully generalized into the entropy Sq = k {I — J2^iPi) / (<? ~ 1) ^ — Sbg)- Mi- 

5^ , croscopic dynamics determines, given classes of initial conditions, the occupation of the accessible 

phase space (or of a symmetry-determined nonzero-measure part of it), which in turn appears to 
determine the entropic form to be used. This occupation might be a uniform one (the usual equal 
CO ^ probability hypothesis of BG statistical mechanics), which corresponds to g = 1; it might be a free- 

scale occupancy, which appears to correspond to g 7^ 1. Since occupancies of phase space more 
complex than these are surely possible in both natural and artificial systems, the task of further 
' fH ' I generalizing the entropy appears as a desirable one, and has in fact been already undertaken in 

^ . the literature. To illustrate the approach, we introduce here a quite general entropy based on a 

^ ' distribution of g- indices thus generalizing Sq. We establish some general mathematical properties 

for the new entropic functional and explore some examples. We also exhibit a procedure for finding, 
given any entropic functional, the g-indices distribution that produces it. Finally, on the road to 
establishing a quite general statistical mechanics, we briefly address possible generalized constraints 
under which the present entropy could be extremized, in order to produce canonical-ensemble-like 



CZ3 ^ stationary-state distributions for Hamiltonian systems. 

. I. INTRODUCTION 

a 

Q ' Since Gibbs' pioneering words [1], it has become well established that standard, Boltzmann-Gibbs (BG) statistical 
mechanics and the associated thermodynamics are valid when certain conditions are satisfied. The typical situation 
occurs for microscopic dynamics exhibiting strong chaos (i.e., positive largest Liapunov exponent) and, consistently, 
^ I exponential mixing, ergodicity, short relaxation times, Euclidean-like occupation of phase space and usual thermo- 
^•s^ I dynamic extensivity. This is the scenario which typically occurs for short-range-interacting many-body Hamiltonian 
0\| ' systems, even if, rigorously speaking, the strict necessary and sufficient conditions remain to be proved. On the 
CO other hand, a vast class of natural and artificial systems exists for which the largest Liapunov exponent vanishes, 
CN situation which is referred to as weak chaos. Weak chaos is typically associated with power-law, instead of exponen- 
tial, sensitivity to the initial conditions and relaxations, (multi)fractal occupation of phase space and thermodynamic 
nonextensivity. 

During the last 15 years a great deal of attention has been focused on systems where the usual BG statistical 
mechanical concepts prove inadequate. Phenomena which exhibit such anomalous behavior occur in systems involv- 
ing long-range interactions (e.g., gravitation), spin-glasses, manganites, long-range micro- or mesoscopic memory, 
turbulence in nonneutral plasma and classical fluids. Levy anomalous diffusion, granular systems, phonon-electron 
'■^ ' anomalous thermalization in ion-bombarded solids, solar neutrinos, cosmic rays, galactic peculiar velocities, cosmo- 
, logical systems, high-energy collisions of particles, black holes, quantum entanglement, lattice Lotka-Volterra growth 
O I model, Hydra viridissima, financial indices and others (see [2, 3, 4] for recent reviews). 

In 1988, a possible generalization of BG statistical mechanics was proposed [5] on the basis of the following entropy 



o 



XI' 1 - 



S =k {qen- Sl = SBG = -ky^P^\np{), (1) 

where fc is a positive constant, usually Boltzmann constant, but from now on taken to be unity for simplicity. This 
generalization of BG statistical mechanics is usually referred to as nonextensive statistical mechanics. Its connection 
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with thermodynamics was estabhshcd in [6] and later redefined in [7] . Its denomination nonextensive comes from the 
following property: if we have two probabilistically independent systems A and B, i.e., Pij{A + B) = Pi{A)pj{B), we 
straightforwardly verify that 

Sg{A + B) = S,{A) + S,{B) + (1 - q)Sg{A)Sg{B) . (2) 

Consequently, q — 1 (the BG case) corresponds to extensivity, whereas q < 1 {q > 1) corresponds to superextensivity 
(subcxtensivity). where the nonnegativity of Sq has been takcni into account. Wc should stress at this point that 
the denomination "nonextensive" can be misleading. Indeed, when special correlations are addressed (instead of 
independent probabilities), it has been shown [8] that Sq is extensive for the appropriate value of q. When no 
scale-invariant (or similar) correlations are present, the appropriate value of q equals unity. 

Within nonextensive statistical mechanics, many of the above cited anomalous systems [2, 3, 4] have found a frame 
of interpretation. The success of such type of approach has led to even more general formalisms. Such is the case of 
Beck-Cohen superstatistics [9, 10, 11, 12], based on fluctuations of parameters such as the temperature [13, 14]. Other 
examples, such as the early proposals by Curado[15] and by Anteneodo and Plastino [16], as well as, more recently, 
Kaniadakis et al formalism based in the K-entropy [17], do exist [18, 19] that have started to exhibit their usefulness 
for applications. 

In this spirit we propose here a self consistent and hopefully applicable generalization of standard non-extensive 
statistical mechanics. We will base our effort on the idea of using not just one q cntropic index but a whole dis- 
tribution of them. The starting point, successfully applied in re-association of folded proteins [20], cosmic rays [21] 
and economics [22], is an ordinary differential equation whose solution exhibits a crossover from a (/-statistics to a 
g'-statistics. This solution might be thought as a generalized "exponential" , whose inverse function is therefore a 
generalized "logarithm" . Two natural manners appear then for extending this type of crossover to a whole spectrum 
of g's, the solution of the full-spectrum ordinary differential equation being interpreted as a quite general exponential 
form, whose inverse is therefore a quite general logarithm. The first manner consists in using the generalized loga- 
rithm for defining a generalized entropy which contains Sq as a particular case; from there, as usual, we obtain the 
stationary-state distribution by optimizing the generalized entropy under appropriate constraints. This distribution 
does not necessarily have the generalized exponential form, a nice property which however does occur for nonextensive 
statistical mechanics. The second manner is to look for a generalized entropy which, when optimized under appropri- 
ate constraints, precisely provides the given generalized exponential. Such entropy will in general not be based on the 
corresponding generalized logarithm, a nice property which, as already said, does occur for nonextensive statistical 
mechanics. The first of these two manners is, although far from trivial, mathematically simpler. This is the one to 
which the present paper is dedicated. The second manner remains as a task to be undertaken in a later occasion. 

In Section II we introduce our formalism by generalizing the logarithmic and exponential functions, as well as the 
corresponding entropic functional. In Section III we study some of their basic properties. In Section IV we illustrate 
the theory with some nontrivial examples, namely a Gaussian distribution of q-exponents and a sum of two delta 
functions. In Section V we introduce and work out the idea of the inverse transformation which enables to find, for 
a given entropic functional, the distribution of exponents that produces it; an illustration is provided as well. In 
Section VI we address the problem of the constraints that can be used in order to obtain a canonical-ensemble-like 
stationary distribution. In Section VII we investigate the relation between the present formalism and Beck-Cohen 
superstatistics. Finally, in Section VIII, we recapitulate our results and discuss some future perspectives. 



II. GENERALIZING THE NONEXTENSIVE ENTROPY FUNCTIONAL 



As well known, BG statistical mechanics is based on the entropy functional 

w 



SBG = J2P^ln^ (3) 



This is a convenient manner of writing Sbg- Indeed, at equiprobability we have pi = l/W, hence the functional takes 
the value In W since the number of (equal) terms precisely cancels Pi in front of the logarithm. The In x function is 
the inverse of the expa; function, which is the solution of the following differential equation [2, 3]: 



y{0) = 1 (4b) 
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In the frame of noncxtensive statistical mechanics Eqs. (3) and (4a) are generalized as follows (Eq. (4b) is 
maintained as it stands): 



w ^ 
^9 = I^Piln<,- 



i=l 



and 



dx 



(5) 



(6a) 



The IngX = [x^~^ — 1]/(1 — q) function is the inverse of the exp^x = [1 + (1 — q)xY/^^~'^^ function, which is the 
solution of Eq. (6a). For g = 1 wc have Ini x = Inx and expj^ x = expx, consequently we fall back to the BG case. 
For g < 1, expq x is taken to be zero if 1 + (1 — q)x < 0. 

Following the above rational used in generalizing BG statistics to nonextensive statistics we can go further and 
examine the effect of the contribution of a whole distribution of q exponents. Within the present generalization the 
entropy is assumed to be of the following form: 



'{/} 



w 1 ^ 



{/} 



(7) 



where the In^jj. x function is the inverse of the exp|yj. x function, defined as the solution of 

dy '■^^ 



dx 



f{r)y^dT 



(8a) 



2/(0) = 1 



with 



/+00 
/(T)dT = l (/(r) > 0, Vr e 7e). 
-OO 



(8b) 
(9) 



We will call the nonncgative, normalized distribution / the q-spectral function ( QSF) since it represents the spectrum 
of q entropic indices that contribute to the entropic functional. In this generalized framework, for /(r) = 5{t — 1) we 
recover the BG entropy Sbg and for /(r) = 5{t — q) we recover the nonextensive entropy Sq, where 5 denotes Dirac's 
distribution. In fact, it is not necessary to impose normalization onto /(r); it is enough to only ask for the integral 



/-oo /(^) finite. If the function /(t) is not normalized, it introduces a very simple modification, namely it 

multiplies the associated generalized logarithm by a constant [23]. It is therefore just for simplicity that we hereafter 
demand /(r) to bo normalized. 

We will now try to establish the form of ln{/} x . By definition of the exp^^j x function, it is 



d exp r f \ X r+oo 

= / m 

J — oo 



dx 



dr 



hence, by setting x = ln{/} y 



hence 



dy 



d[ln{f}y] 



/+00 
f{r)y^dT, 
-oo 



X 



At equiprobability (i.e., pi = 1/W^) we have 



f(T)u'dT\ du (Vx e (0,+oo)). 



(10) 



(11) 



(12) 



'{/} 



-w 



/(r) u^dT 



du 



(13) 
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III. PROPERTIES OF THE GENERALIZED LOGARITHM AND EXPONENTIAL FUNCTIONS AND 

THE CORRESPONDING ENTROPY 

Wc list now some useful properties that can be easily proven for any normalized distribution f{q), given the fact 
that f{q) is nonnegative. 

ln{/} 1=0, (14) 

hence 

exp{/} = 1. (15) 

Also 



^ 1 



^exp^/}^ 



= 1 , (16) 

x=0 



as well as monotonicity, more precisely 

dx 

and 



ln{/} X > 0, Va; e (0, +oo) , (17) 



^ ®^P{/} X>0,WXG Aexpyy , (18) 

where <E Aexpif^ is the set of admissible values of x for the nonnegative expjj-j. a; function. When no negative q 
contributes (i.e., if f{q) = 0, ^q < 0), then the following properties hold also: 

^— 2 X < {concavity) , (19) 



and 

dx^ 



exp|y} a; < (convexity) . (20) 



Analogously, when no positive q contributes (i.e., if f{q) = 0, Vg > 0), then 

dx'^ 



d? 

— 2 X > {convexity) , (21) 



and 

e^P{f} X > (concavity) . (22) 
Also, if no q above unity contributes (i.e., if f{q) = 0, Vg > 1), then 

lim In r f 1 X — +00 , (23) 

x—>+oo ^ ' 

and, if no q below unity contributes (i.e., if f{q) = 0, < 1), then 

li™ ^^{f} x = -oo. (24) 

X— >0+ 
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IV. EXAMPLES OF THE GENERALIZED ENTROPY 
A. Gaussian distribution 

As a first example of QSF, let us consider a Gaussian distribution of entropic indices around a central index q with 
a variance cr > 0: 

1 jT-lf 

/(r) = ^=e-^ (25) 
Using equation (12) we obtain, for the corresponding generalized logarithm: 



f I ^ ~ «^ /I - g - CT^lna; 

erf -^=- ] - erf 



\Plu ) \ ^/2a 



(26) 



Where erf is the error function. As special cases we obtain: 



— 1 

^Hifi}^= ' (27) 



and 



The entropy functional is given by 



' ,/l-g + aMn(l/ft) 



\/2(7 



(29) 



In Fig. 1 we present, for typical values of (q, a), the entropic function S{<j,cr}(p), the entropy S{q^„}{p) for a two-state 
system, and the microcanonical entropy S^q^^^iyV) as a function of the number of states W. 

B. Two entropic indices 
The general problem of two indices in a QSF corresponds to the form: 

/(r) = ai^(r - gi) + a25(r - q^) (30) 

(If we wish to have a normalized QSF wc need to require 01+02 = 1.) Since the general case of Eq. (30) is 
quite untractable, we will focus our attention onto two special cases, referred to as the symmetric and the asymmetric 
ones. Let us first consider the case where the two indices are symmetric with regard to g = 1, i.e., qi = 1 — L/2 and 
92 = 1 + L/2 and a\= a2 = 1/2. We then have: 

f{r) = \ [5{t - 1 + L/2) + 5{t-1- L/2)] (31) 

hence, using (12), 

r 2 , Aarctan{x^/'^)--K 

In the limit L — > 0+ we recover Inj;, a; = In x and we will thus be driven back to the BG case. The entropy functional 
is given by: 
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w 



Aarctanip^ ^^^) — tt 



(33) 



It can be readily seen that, for equiprobability, the entropic functional is given by 

4arctan{W'^^^) - tt 
= ;^ 



(34) 



Notice that, when W oo, the entropy Sl remains finite, in contrast with the BG case. 

In Fig. 2 we present, for typical values of L, the entropic function sl{p), the entropy Sl{p) for a two-state system, 
and the microcanonical entropy Sl{W) as a function of the number of states W. 

We consider now the asymmetric case, defined as the special case of Eq. (30) where one of the indices equals unity. 
Then we have: 



/(t) = aiS{T - 1) + aqS{T - q) (ai + a, = 1) 



(35) 



This form might be of physical relevance. As already mentioned, unless ai + Uq = 1, /(r) is not normalized (see 
[23]). The differential equation corresponding to the (ai, a^) generic case has already been successfully used (in order 
to describe, for instance, the re-association of folded heme proteins [20] and the flux of cosmic rays [21]). 

It can be proven that the corresponding generalized logarithm is: 



In 



ai(l - q) 



In 



ai 



ai + aq 



(l-x-«) 



(36) 



We straightforwardly verify that iT^{ai,o,q} x = \nx and that ln{o,i,g} x = In^ x. The corresponding entropic functional 
can be writen as: 



'{ai,a<,,g} 



W 

In the microcanonical case this can be writen as: 



^ W r 

— ^^pan i-^(i-pr) 

ai(l - g) ai+ aq \ J 



(37) 



ai(l-g) 



In 



1 - 



ai 



ai + aq 



(1 - w^-'i) 



(38) 



Like in the previous case, it can be shown that, for q> \. the entropy remains finite when W +qo. In Fig. 3 we 

present, for typical values of {q,aq) and Oi = 1 — a,, the entropic function .S{ai,a,.g}(p); the entropy S^a^ a,j,q}{p) 
two states, and the entropy S^ai,ag,q}{W) in the microcanonical ensemble as a function of the number of states W. 

V. CALCULATION OF THE QSF FOR A GIVEN ENTROPIC FUNCTIONAL 

We have seen so far how from a given QSF we can produce the corresponding entropic functional. We will now work 
in the reverse way: given a specific entropic functional, we will find (if possible) the QSF that produces it. Consider 
a general entropic functional of the form: 



w 



(39a) 



s{x) = X ln{/} - 



(39b) 
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We have: 

, . du 

hirnx = 



{f}^= / 7+^ 



1 /_oo f{T)u^dT 



— (mrfi x) = ,^ 



/ 
/. 



+ 00 

f{T)x-dT = ^— ^ 



+ 00 1 

r Inx J_ 



We set: 



Then, from equation (40), we obtain: 



/(r)e- "^-dr = ^ (40) 



a; = — i Ina; (41) 



1 r + oc 1 

^ / f{T)e"^^dT = ^= ■ -r-, ^ (42) 

The LHS of equation (42) is however nothing but the Fourier transform of /. Thus, inverting the transform we 
have: 

fil) = ir / ^TT ^ ^ (43) 

Inserting the entropy functional of Eq. (7) into equation (43) wc finally get: 

^^'^ = ^L .(e--)-.^Ke-)]'^ ^''^ 

Equation (44) is quite important. Indeed, it shows that the QSF corresponding to a large class of entropy functionals 
can be explicitly calculated. It is straightforward to check that for s given by BG or by nonextensive statistics we get 
f{q) = S{q — qo) as anticipated. 

In order to be able to derive the normalized QSF associated with a given entropy, the entropy functional must fulfill 

the following requirements. 

1. It must be possible to write the total entropy 5 as a sum of the entropic function s for each state (Eq. (39a)). 

2. The function s must satisfy s(l) = 0, which is in fact a quite reasonable requirement for an entropy. 

3. Furthermore, we must have: 



ds{x) 



dx 



= -1 (45) 



This condition is equivalent to having the QSF normalized to unity. If we abandon the normalization of the 
QSF then we can consistently drop this last requirement. 



8 



4. The function s{x) must be defined (or analytically continued) on the unitary circle and it must also be differen- 

tiablc in the same domain. 

5. The integral of Eq. (44) must converge. 

As a nontrivial illustration, we will now use the present method to find the QSF associated with an exponential 
entropic form. Let us assume 



i=l 

hence 



w _^ 

^ = l]ft(l-eV), (46) 



s{x) =x[l-e—j . (47) 

It is trivial to see that Eq.(47) fulfills all the criteria set above, and we can thus find a normalized QSF for it. Using 
Eq.(44) we get: 



n=0 

Although different, entropy (46) has some resemblance with that introduced by Curado [15]. We claim no particular 
physical justification for the form (46). In the present context, it has been chosen uniquely with the purpose of 
illustrating the mathematical procedure involved in the inverse QSF problem. 



VI. THE CANONICAL ENSEMBLE 



To obtain, for the canonical ensemble within the standard nonextensive thermostatistics, the distribution corre- 
sponding to the physically relevant stationary state, we must optimize the entropic functional 



w 



Besides normalization of the probabilities, we use the following energy constraint (see [7, 24] and references therein): 



= U, (50) 



where {Ei} are the eigenvalues of the Hamiltonian. The optimizing distribution can be written as a generalized 
exponential form as follows: 

where /3 is an inverse-temperature like parameter. As we see, the stationary distribution involves precisely the inverse 
function of the generalized logarithm that entered into the definition of the entropy. The following question arises: 

In the general case considered in the present paper, is it possible to write the energy constraint in such a way that 
the stationary state distribution is described by the inverse function of the generalized logarithm used in the definition 
of the entropy? Let us right away say that, if this is possible, we have not sueccded in doing it. It seems that this 
nice property is not very general, and it appears to be lost when we extend further (or differently) than nonextensive 
thermostatistics. 
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Let us present an attempt along that line. We assume that the entropy is given by Eq. (7). The corresponding 
energy constraint is assumed to be given as follows: 

In the case where f{q) = S{q — qo) we get back to the usual nonextcnsivc constraint. 

The stationary state canonical probability distribution is the one that optimizes the entropy functional under this 
energy constraint. In all the examples that we have given in Section IV, this calculation is analytically untractable. 
Nevertheless, it can be trivially shown that the solution is not the inverse function of the generalized logarithm that 
entered into the definition of the entropy. It seems thus, that the exponential form is not a generic feature but rather 
a peculiar "coincidence", only valid for the nonextensive family of entropies Sq. It remains, however, as an open and 
rathca- challenging problem to come up with an energy constraint that would produce the corresponding generalized 
exponential in the generic case. 



(52) 



VII. CONNECTION WITH THE BECK-COHEN SUPERSTATISTICS 



Recently, a generalization of the BG statistics has been proposed by Beck and Cohen. This generalized superstatis- 
tics assumes that there may be strong deviations from the classical BG case when there are temperature fluctuations 
within the system. Along that line, the Boltzmann factor ex.p{—pE) is generalized as follows: 



r+oo 

B{E) = / g{(i)e-l'^dp (53) 

Where g{l3) is the probability distribution of the inverse tempature in the system. Quite recently, an entropic 
functional has been derived that corresponds to superstatistics. This functional is of the form: 



Jo 1-^ 



M 



(54a) 



(54b) 



B{y) 

!rB{u)du 



where K{y) = ^ ^ (54c) 



In Eqs. (54), E* stands for the lowest admissible energy for the system and a is a Lagrange multiplier. Using Eq. 

(54b) and Eq. (44) we can in principle get the QSF f{q) for a given temperature distribution fimction g{l3). Thus, in 
principle at least, the various superstatistics can be accomodated into the presently proposed formalism. 

There are certain cases however, where the present formalism can go further. For example, it appears that through 
superstatistics we can only produce the so-called nonextensive entropies Sq iov q>l, while in the present formalism 
we can have them for arbitrary values of q. 



VIII. CONCLUSIONS 



In the present work, we study the properties of a generalized entropic function based on a distribution of q entropic 
indices. We introduce the idea of g-spectral function(QSF) and exhibit three nontrivial examples: the case of a 
Gaussian, as well as both symmetric and nonsymmctric choices within a sum of two delta functions. Wc also introduce 
the notion of an inverse QSF problem, in the sense that, given some entropic functional, we present a procedure for 
explicitly calculating the QSF that produces it. 

Also, we have briefly discussed the canonical ensemble and its corresponding energy constraint. Wc have come to 
the rather unexpected conclusion that, in the general case, the stationary state probability distribution seems to be 
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not the generalized exponential which is the inverse function of the generalized logarithm used to define the entropy. 
Altering the form of the constraint in order to produce such exponential remains therefore as an open problem. 

We have also pointed out the connection of the present work to the Beck-Cohen superstatistics. Incidentally, we 
have shown that there are examples of entropies that can be produced within the current theory which appear to be 
not admissible within the frame of Beck-Cohen superstatistics. 

It would certainly be very interesting to apply the present formalism to concrete physical systems (e.g., high 
precision experiments in turbulent systems) that would defy more conventional approaches. On a different vein, it 
would be theoretically valuable to find nontrivial special cases for which the canonical ensemble distribution could be 
explicitly calculated, since this would allow a deeper analysis of the properties of the optimizing distribution. 
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FIG. 1: The entropic function siq^^y{p) {top), the entropy S^q^a^y{p) for a two-state system {middle), and the VF-state 
microcanonical-ensemble entropy S^q^a}{W) {bottom). The four Unes correspond to g = 1 and a = 0.1, a = 0.5, cr = 1 
and a = (BG case). 




FIG. 2: The entropic function sl{p) (top), the entropy Sl{p) for a two-state system (middle), and the VK-state microcanonical- 
ensemble entropy Sl{W) {bottom). The four lines correspond to L = 1 , L = 5, L = 10 and for L = (BG case). 
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FIG. 3: The entropic function S{ai,a,,q} (p) {top), the entropy S^ai,ag,q}ip) for a two-state system {middle), and the W^-state 
microcanonical-ensemble entropy S^^ai,ag,q}{W) {bottom). The four hnes correspond to g = 2, oi = 1 — Oq, and Uq = 0.9, 
ttq = 0.6, Oq = 0.3 and Oq = (BG case). 



